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Abstract 



The nature of the fuel that drives today's cosmic acceleration is an open and 
tantalizing mystery. The braneworld theory of Dvali, Gabadadze, and Porrati 
(DGP) provides a context where late-time acceleration is driven not by some 
energy-momentum component (dark energy), but rather is the manifestation 
of the excruciatingly slow leakage of gravity off our four-dimensional world into 
an extra dimension. At the same time, DGP gravity alters the gravitational 
force law in a specific and dramatic way at cosmologically accessible scales. We 
derive the DGP gravitational force law in a cosmological setting for spherical 
perturbations at subhorizon scales and compute the growth of large-scale 
structure. We find that a residual repulsive force at large distances gives rise to 
a suppression of the growth of density and velocity perturbations. Explaining 
the cosmic acceleration in this framework leads to a present day fluctuation 
power spectrum normalization as < 0.8 at about the two-sigma level, in 
contrast with observations. We discuss further theoretical work necessary to 
go beyond our approximations to conflrm these results. 



*E-mail: 
^E-mail: 
^E-mail: 



lue@cern.ch 

rsl23@nyu.edu 

glenn . st arkman@cer n . ch 



I. INTRODUCTION 



The discovery of a contemporary cosmic acceleration [1,2] is one of the most profound 
scientific observations of the 20th century. We are now challenged to answer the open and 
tantalizing question of what drives that acceleration. While a conventional explanation 
exists (i.e. dark energy - some new negative- pressure energy-momentum component), an 
intriguing line of thought is gaining attention: the accelerated expansion is not a result of 
yet another ingredient in our already gunky cosmic gas tank, but rather is a signal of our 
lack of understanding of gravitational physics on large scales [3-14]. 

Being able to observationally differentiate the two possibilities, dark energy versus mod- 
ified gravity, is an essential component in developing the modified-gravity paradigm. One 
can easily envision some modified-gravity model leading to an expansion history that can 
be identically reproduced by some dark-energy model. Thus, observations that depend only 
on anomalous expansion histories are insufficient to tease out the acceleration's root cause. 
However, as we argue in [15], if one attempts to modify cosmology at today's Hubble scale, 
Hq, through altering the equations governing gravitational dynamics, then generically one 
expects that the gravitational force law of an isolated mass source is altered even at distance 
scales much smaller than Hq^. This effect can then be exploited to differentiate between 
a modified-gravity explanation of today's cosmic acceleration and dark energy (where the 
gravitational force law remains unaltered) [16,15]. 

A leading contender in modified-gravity explanations of acceleration is the braneworld 
theory of Dvafi, Gabadadze, and Porrati (DGP). In this theory, gravity appears four- 
dimensional at short distances but is altered at distances large compared to some freely 
adustable crossover scale Tq through the slow evaporation of the graviton off our four- 
dimensional braneworld universe into an unseen, yet large, fifth dimension [17-19]. DGP 
gravity provides an alternative explanation for today's cosmic acceleration [3,4]: just as 
gravity is conventional four-dimensional gravity at short scales (r <^ Tq) and appears five- 
dimensional at large distance scales (r » tq), so too the Hubble scale, H{t), evolves by the 
conventional Priedmann equation at high Hubble scales but saturates at a fixed value as 
H{t) approaches r^^. Thus, if one were to set ro ~ -f^o"^, where Hq is today's Hubble scale, 
then DGP gravity could account for today's cosmic acceleration in terms of the existence of 
extra dimensions and a modification of the laws of gravity. The resulting cosmic expansion 
history is specific and may be tested using a variety of cosmological observations [4,20-22]. 
However, can we distinguish between DGP gravity and a dark energy model that mimics 
the same cosmic expansion history? 

We would naively expect not to be able to probe the extra dimension at distances much 
smaller than the crossover scale ro = Hq^. However, in DGP, although gravity is four- 
dimensional at distances shorter than tq, it is not four-dimensional Einstein gravity - it is 
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augmented by the presence of an ultra-light gravitational scalar. One only recovers Einstein 
gravity in a subtle fashion [23-26], and a marked departure from Einstein gravity persists 
down to distances much shorter than tq. For example, for tq ~ H^^ and a central mass 
source of Schwarzschild radius r^, significant and cosmologically-sensitive deviations from 
Einstein gravity occur at distances greater than [25-28] 



Thus a marked departure from conventional physics persists down to scales much smaller 
than the distance at which the extra dimension is naively hidden, or for our discussion here, 
the distance at which the Friedmann equation was modified to account for accelerated cosmic 
expansion. This alteration of gravitational interactions provides a way of differentiating 
between DGP gravity and dark energy models, and is consist with the argument we put 
forth in [15]. 

Imminent solar system tests have been shown to be capable of probing the residual devia- 
tion from four-dimensional Einstein gravity at distances well below r* [27,28]. Nevertheless, 
it would be ideal to test gravitational physics where dramatic differences from Einstein grav- 
ity arc anticipated. A detailed study of large scale structure in the Universe can provide such 
tests of gravitational physics at large distance scales. Unfortunately, prior analyses related 
to modified-gravity explanations of cosmic acceleration [16,15] are not applicable here. The 
modified force law is, in effect, sensitive to the background cosmological expansion, since this 
expansion is intimately tied to the extrinsic curvature of the brane [3,29], and this curvature 
controls the effective Newtonian potential. A more careful analysis must be performed. In 
the next section we briefly review DGP gravity and identify the new force law necessary 
to calculate how large scale structure evolves in this cosmological model. We then proceed 
and compare those results to the standard cosmology, as well as to a cosmology that exactly 
mimics the DGP expansion history using dark energy. Finally, we discuss the observational 
implications of our results on the growth of structure in DGP gravity and conclude with 
some remarks and a discussion of future work needed to improve upon our treatment. 



We review the important points of the DGP braneworld model, including cosmology. 
We then modify the calculation performed in [27] to do determine the gravitational force 
law in an evolving cosmological background, rather than in a static background deSitter 
space. This calculation is the new result in this paper that allows one to answer questions 
of cosmological interest developed in the next sections. 




(1.1) 



II. DGP GRAVITY 
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A. The DGP Model 



Consider a braneworld theory of gravity ( one in which ordinary particles and fields, other 
than the graviton, are confined to a three-dimensional hypersurface - the brane - embedded 
in a higher dimensional space - the bulk) with an infinite-volume bulk and a metastable 
brane graviton [17]. We take a four-dimensional braneworld embedded in a five-dimensional 
Minkowski spacetime. The bulk is empty; all energy-momentum is isolated on the brane. 
The action is 



J d'xy/^ R + J d'xyf^) Cm + Sgh . (2.1) 



IGtt 

The quantity M is the fundamental five-dimensional Planck scale. The first term in Eq. (2.1) 
corresponds to the Einstein-Hilbert action in five dimensions for a five-dimensional metric 
Oab (bulk metric) with Ricci scalar R. The term Sgh is the Gibbons-Hawking action. In 
addition, we consider an intrinsic curvature term which is generally induced by radiative 
corrections by the matter density on the brane [17]: 

-^M| / d'x^f^) . (2.2) 

Here, Mp is the observed four-dimensional Planck scale (see [17-19] for details). Similarly, 
Eq. (2.2) is the Einstein-Hilbert action for the induced metric g^J on the brane, R^*"^ being 
its scalar curvature. The induced metric is^ 

= d.X^d^X^g^B , (2.3) 

where X^{x^^) represents the coordinates of an event on the brane labeled by x^. The action 
given by Eqs. (2.1) and (2.2) leads to the following equations of motion 

:^Gab + S{hrane)G^^l = ^T^slbrane , (2.4) 
zr iVip 

where Gab is the bulk Einstein tensor, is the Einstein tensor of the induced brane 
metric, and T^b| brane is the matter energy-momentum tensor on the brane. We have defined 
a crossover scale 

-■0 = ^ ■ (2-5) 

This scale characterizes that distance over which metric fluctuations propagating on the 
brane dissipate into the bulk [17]. 



^Throughout this paper, we use A,B,... = {0, 1, 2, 3, 5} as bulk indices, /x, z/, . . . = {0, 1, 2, 3} as 
brane spacetime indices, and = {1, 2, 3} as brane spatial indices. 
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B. The Cosmological Background 



Let us review some important details of the cosmological background, for a general ex- 
pansion of a spatially-flat brane into a Minkowski-fiat bulk. We are primarily interested 
in the late evolution of the Universe, in particular the matter-dominated era where the 
energy-momentum content of the universe is wcll-rcprcscntcd by a prcssurclcss distribu- 
tion of galaxies, spatially homogeneous on the largest scales. The spatially homogeneous 
cosmological background of such a Universe is driven by energy-momentum given by 

rilbrane = S{z) diag (p(t), 0, 0, 0, 0) , (2.6) 

with spacetime geometry dictated by the line element [3] 

ds^ = (l T dr^ - (l T ^\z\^ [SijdX'dX^] - dz" . (2.7) 

Here, dot refers to differentiation with respect to the cosmological time r, the coordinate 
z is the extra dimension, and the brane scale factor, a(r) satisfies a modified Friedmann 
equation 

where H{t) — a/d. The two choices of sign represent two distinct cosmological phases. The 
phase of interest (the self-accelerating phase) corresponds to the lower sign, but we keep 
both for the sake of completeness. 

This new Friedmann equation Eq. (2.8) already makes the theory distinct from standard 
ACDM cosmology, and observational signatures constraining DGP cosmology have been 
considered in e.g. [4,21,22]. Using constraints from Type lA supernovae [21], the best fit Tq 
is ro = 1.2llg:^^ifo~\ where Hq is today's Hubble scale. Taking Hq ^ 70 km s~^Mpc~\ it 
implies ro ~ 5 Gpc. 

However, we wish to focus on those properties of DGP gravity that are affected by the 
modification of the force law, and cannot be mimicked by some dark energy component 
that produced the same expansion history. Following this program, we focus particularly at 
distance scales much smaller than the Hubble radius, H~^. As described in the introduction, 
the gravitational force law is significantly different from four-dimensional Einstein, even 
at these short distance scales. We wish to determine the form of the corrections in the 
background of the expected matter-dominated cosmology, Eq. (2.8). 

We are concerned with processes at distances, r, such that rH <^ 1. Under that circum- 
stance it is useful to change coordinates to a frame that surrenders explicit brane spatial 
homogeneity but preserves isotropy 
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r(T, X') = a(r)A 



for all z and where = SijX^X^. The hne element becomes 



1t2{H + H/H)\z\ - {H' + Hy dr - [l ^ 2H\z\] (1 + H'r')dr' + r'dil 



(2.9) 
(2.10) 

— dz"^ , 
(2.11) 



where here dot repreresents differentiation with respect to the new time coordinate, t. More- 
over, H — H{t) in this coordinate system. All terms of 0{r^H^) or 0{z'^H'^, zHrH) and 
higher have been neglected. 

The coordinate system in Eq. (2.11) will be the most useful when considering the cosmo- 
logical scenarios we are interested in. One can see that the bulk is a Rindler space. This has 
a fairly natural interpretation if one imagines the bulk picture [3,29]. One imagines riding a 
local patch of the brane, which appears as hyperspherical surface expanding into (or away 
from) a five-dimensional Minkowski bulk. This surface either accelerates or decelerates in 
its motion with respect to the bulk, creating a Rindler-type potential. 



C. Nonrelativistic Matter Sources 



We are interested in finding the metric for compact, spherically symmetric overdensities 
in the background of a matter-dominated cosmology. Because we are only concerned with 
distance scales such that rH <g 1, then to leading-order in r^i?^ and zH, the solutions to 
the field equations Eqs. (2.4) are also solutions to the static equations, i.e. the metric is 
quasistatic, where the only time dependence comes from the slow evolution of the extrinsic 
curvature of the brane. To be explicit, we are looking at the nonrelativistic limit, where the 
metric, or the gravitational potentials, of a matter source depends only on the instantaneous 
location of its elements, and not on the motion of those elements. 

Under this circumstance, one can choose a coordinate system in which the cosmological 
metric respects the spherical symmetry of the matter source. Let the line element be 

ds^ = N\t, r, z)dt^ - A^{t, r, z)dr^ - B^t, r, z)[de^ + sin^ Odcf)'^] - dz^ . (2.12) 

We are interested in small deviations of the metric from flatness so we define new functions 
{n{t, r, z), a{t, r, z), b{t, r, z)} such that 

N{t,r,z) = 1 +n{t,r,z) 

A{t,r,z) = l + a{t,r,z) (2.13) 
B{t,r,z) = r[l + b{t,r,z)] . 
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The key is that because we are interested primarily in phenomena whose size is much smaller 
than the cosmic horizon, the effect of cosmology is almost exclusively to control the extrinsic 
curvature, of the brane. This can be interpreted as a modulation of the brane 's stiffness or 
the strength of the scalar gravitational mode. 

We take the energy-momentum tensor to now be 

r^lbrane = S{z) diag + r) , 0, 0, 0, 0) , (2.14) 

where the source mass is an overdensity of compact support (i.e. its extent is some ra- 
dius, R <^ H~^). Given a source mass whose overdensity with respect to the cosmological 
background goes as 5p{r,t), we may define an effective Schwarzschild radius 

^9ir,t)^j^l r'Sp{r,t)dr. (2.15) 

We solve the perturbed Einstein equations in quasistatic approximation by generalizing the 
method used in [27], obtaining the metric of a spherical mass overdensity Sp{t,r) in the 
background of the cosmology described by Eq. (2.11) (rather than dcSittcr space). The 
metric on the brane, using the residual gauge fixing b{t,r)\z=o = 0, is then given by 

rn'it, r)U=o = ^ [1 + ^(0] - + (2.16) 
a(t,r)|,=o = §[l-A(r)] + ^i/V (2.17) 

where dot now denotes differentiation with respect to t and prime denotes differentiation with 
respect to r. Note that the background contribution is included in these metric components. 
The quantity A(r) is defined as 



A(rJ 



(2.18) 



and 



Though it is not of explicit interest here, the full 2;-dependence of the metric may be deduced 
from Eqs. (2.16) and (2.17) using equations laid out in the appendix of Ref. [27] with trivial 
alterations accounting for the differing cosmological background. 

The result Eqs. (2.16)-(2.19) is valid for r <^ and r ^ tq, but only if the spatial sup- 
port of Sp{r, t) extends only to radii much less than H~^, so that there is a clear distinction 
between the matter making up the overdensity and the cosmological background. The result 
is virtually identical to the strictly static deSitter background case except there /? = l±2roH. 
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One may also confirm that in the absence of perturbations (i.e., 5p or i?g = 0), the back- 
ground metric Eq. (2.11) is a consistent, quasistatic solution. This point is analogous to 
the well-known idea that one may reproduce the Priedmann equation in matter-dominated 
cosmologies with just the Newtonian interaction bewteen matter particles. 

From inspection of Eq. (2.19), we see that, in addition to Tq, there exists a new transition 
scale 

1/3 

(2.20) 



such that when r <^ r*, the Einstein phase, the metric functions on the brane reduce to 



n — — 




When r r^, the weak- brane phase, the metric functions on the brane become 



n — — 



Rg^ 

2r 



1 + 



1 

3^ 



2r 



1 



2^ 

2 



(2.21) 
(2.22) 

(2.23) 
(2.24) 



In this phase, the extra scalar mode, the would-be radion, alters the effective Newton's con- 
stants for the gravitational potentials represented by n{t, r) —n{t, r)| background) the Newtonian 
potential, and a{t,r) — a(i, r)|background, the gravitomagnetic potential. 

One may simply check that the full {t, r, 2;)-dependent metric satisfies the complete mod- 
ified Einstein equations Eqs. (2.4) to the desired order. Note that to this order of precision, 
the velocity of the matter distribution Sp{t,r) docs not affect spacetime geometry (until 
order v'^ or vrH) , thus corroborating the quasistatic approximation. 



D. Caveats 



The approximations v <^ 1 {v represents peculiar matter velocities) and rH ^ 1 play 
a role in several places and allow a scries of crucial simplifications that need to be spelled 
out. These two approximations are lumped together because the Hubble-flow velocity and 
peculiar velocities play almost identical roles in the relevant field equations. The following 
are the operational simplifications: 

• Nonrelativistic, quasi-static sources. Source-velocity dependent contributions to the 
gravitational field are subleading. One may use the static Einstein equations and still 
be assured that the metric on the brane is accurate to 0{v'^ ,r'^H'^). 
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• Near-field regime. Related to the above simplification, the source evolves slowly enough 
that radiative effects are negligible at these radii. We may safely avoid scalar radiation 
on the brane and gravity-wave evaporation into the bulk and other radiative bulk 
effects [30,29]. These radiative scalar modes may also have classical instabihties that 
might become relevant outside this regime [31]. 

• Resolving background matter from perturbation. This simplification is specific to 
DGP gravity. The metric components Eqs. (2.16) and (2.17) depends on the back- 
ground Hubble expansion, and the calculation crucially depends on the assumption 
that the overdensity does not alter the background cosmology. This can been seen to 
be self-consistent in the DGP field equations only when the radius of support for the 
overdensity is much smaller than the Hubble radius, H~^. 

• Geodesic motion is Newtonian, i.e. the geodesic equation reduces to Newton's second 
law where the potential is goo or n{t, r). 

If we stray too far from the assumptions f <^ 1 and rH ^ 1, then effects safely disre- 
garded may start intruding into and complicating the analysis, particularly when H Ho, 
introducing additional effects of equal significance to the ones included here. 



Let us consider the evolution of a spherical top-hat perturbation 5{t, r) of top-hat radius 
Rt, where p{t,r) — p{t){l -\- 5) is the full density distribution and p{t) is the background 
density. At subhorizon scales {Hr <S 1), the contribution from the Newtonian potential, 
n{t,r), dominates the geodesic evolution of the overdensity. Prom Eq. (2.23) it follows the 
equation of motion for the perturbation 5 is. 



where the definition of e = Sr^Rg/d^'^Rf follows from the identification of the expression 
in square brackets with 1 + A(r) (see Eq. (2.18)), and we have restricted ourselves to the 
self- accelerating branch (i.e., the lower sign choice in all equations in the previous section). 

For clarity, we may recast the time evolution of (3 and e in terms of S and the time- 
dependent value of flm- Defining 0,^{t) = n^Vff2u) , and using the Friedmann equation. 



III. GROWTH OF DENSITY PERTURBATIONS 



A. Spherical Perturbations 




(3.1) 



Eq. (2.8), 
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We stress that Jl^ is a time-dependent quantity - it goes to unity at high redshift 1 <^ z <^ 
Zeq, where the evolution is Einstein-deSitter (but in the matter dominated regime), and at 
present it reduces to the usual value that we denote as O,^ — flmi^ — 0). We see that P is 
negative, of order unity at present, and approaches minus infinity at high redshift, whereas 
e is proportional to with a coefficient of order unity. Note that for large 6, Eq. (3.1) 
reduces to the standard evolution of spherical perturbations in general relativity. However, 
when S is small, the correction term in the square brackets can be noticeably different from 
unity. 

B. Linear Growth 

Let us focus first on linear perturbation growth at scales r <^ H~^. In this regime 
S{r,t) <^ 1, therefore one is always in the weak-brane regime,^ r 3> r^,, and the only effect 
of DGP gravity is a modification of Newton's constant. Equation (3.1) reduces to 

5 + 2HS = AnGp (^1 + ^) ^ ■ (3-3) 

Note that the effective Newton's constant, 

G^ = g{i^1}j, (3.4) 

is time-dependent. Since /3 is negative, as time goes on the effective gravitational constant 
decreases, and this extra repulsion (compared to general relativity) leads to suppressed 
growth. For example, ff = 0.3, Geff/G = 0.72, 0.86, 0.92 at ^ = 0, 1, 2. 

The growing-mode solution of Eq. (3.3), D+, is shown as a function of redshift z in Fig. 1. 
The top panel shows as dashed fines the ratio of in DGP gravity to that in a dark energy 
(DE) scenario with the same Friedmann equation but standard gravity, for two values of the 
present matter density = 0.3 (top) and — 0.2 (bottom). Notice how the change in 
the effective Newton constant leads to a suppression of Incidentally, this suppression 
is about two times larger than for models of modified gravity (with the same expansion 



^This is always true for top-hat perturbations, but in practice the size r is related to the amplitude 
6 through the perturbation spectrum. However, for la fluctuations of scale r = 10 — 100 Mpc/h 
with typical profiles given by the two-point correlation function, r* corresponds to 5 — 15 Mpc/h. 
Therefore perturbations accessible to large-scale structure surveys are a natural probe of DGP 
gravity in the weak-brane regime. 
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FIG. 1. The top panel shows the ratio of the growth factors (dashed Unes) in DGP gravity 
[Eq. (3.3)] and a model of dark energy (DE) with an equation of state such that it gives rise 
to the same expansion history (i.e. given by Eq. (2.8), but where the force law is still given by 
general relativity). The upper line corresponds to = 0.3, the lower one to 0^ = 0.2. The solid 
lines show the analogous result for velocity perturbations factors /. The bottom panel shows the 
growth factors as a function of redshift for models with different expansion histories, corresponding 
to (from top to bottom) ACDM {Vt%^ = 0.3), and DGP gravity with Vt^^ = 0.3,0.2 respectively. 

history) that obey the Birkhoff's law [16,15]. The lower panel compares the growth factor 
Dj^ to that in the standard cosmological constant scenario (with Vt^^ = 0.3 and Q\ = 0.7), 
again for Q'^ = 0.3 (top) and Q'^ = 0.2 (bottom). We see here that the change in the 
expansion history (from a cosmological constant to DGP) leads to an additional suppression 
of the growth. In the language of dark energy, this is because the non-standard term in 
the Friedmann equation Eq. (2.8) can be thought of as a contribution from a dark energy 
component with an effective equation of state given by 

^^eff = - , (3.5) 

i T i I'm 

therefore, for fixed such a term dominates the expansion of the universe earlier in DGP 
gravity than in DE models with a cosmological constant, leading to an enhanced expansion 
rate H and therefore and additional suppression over the one provided by the change in the 
force law. We will examine the observational consequences of this in Sec. IV. 

The growth of velocity perturbations is also a useful observable, and it follows directly 
from the continuity equation in the linear approximation. It is specified by / = In /din a 
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FIG. 2. Numerical solution of the spherical collapse. The left panel shows the evolution for a 
spherical perturbation with 5j = 3 x 10~^ at Zi = 1000 for = 0.3 in DGP gravity and in ACDM. 
The right panel shows the ratio of the solutions once they are both expressed as a function of their 
linear density contrasts. 



and it can be parametrized in terms of the time variable VL^ 
differential equation for f{flm) directly from Eq. (3.3), 



in fact, one finds the following 



df 



dVtr 



+ 



1 



3fi™('i - n 

r 



(2-^]^)/ + (i + a^)/^ 



l(l + r]^)(l + 21]; 

3 + 



(3.6) 



whose numerical solution follows approximately f{i^rn) — ^m^' which can be contrasted with 
the standard /(^m) — models with a cosmological constant. The top panel 

of Fig. 1 shows the ratio of / for DGP and DE models with the same expansion history, 
for = 0.3, 0.2, showing that differences of at least 10% are expected, whereas comparing 
DGP to cosmological constant models with the same larger differences are obtained, e.g. 
for = 0.2, jDGP^jACDM _ Q_g3^ These deviations are well within the range that can be 
probed with current redshift surveys. 

C. Non-Linear Grow^th 



The left panel in Fig. 2 shows the full solution of Eq. (3.1) with an initial condition of 
Si = 3 X 10^^ at Zi = 1000 for = 0.3, and the corresponding solution in the cosmological 
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constant case. Whereas such a perturbation collapses in the ACDM case ai z — 0.66 when its 
hnearly extrapolated density contrast is 6c — 1.689, for the DGP case the collapse happens 
much later at z — 0.35 when its 5c — 1.656. In terms of the linearly extrapolated density 
contrasts things do not look very different, in fact, when the full solutions are expressed as 
a function of the linearly extrapolated density contrasts, Sn^ — D+5i/{D+)i they are very 
similar to within a few percent (right panel in Fig. 2) . This implies that all the higher-order 
moments of the density field will be very close to that in ACDM models. Indeed, such 
moments are determined by the vertices i/n defined from {i/i = 1) 

S{Si^) = E ^ ^L, (3.7) 

n=l ^■ 

e.g. the skewness is = 3^2 [32], up to smoothing corrections that depend on the trans- 
formation from Lagrangian to Eulerian space. However since = 31^5 /dSy^^ at 5\i^ = 0, the 
f„'s in DGP gravity will all be very similar to those in ACDM (we have checked this explic- 
itly for 5*3, obtaining less than 1% change). This can be useful in the sense that it allows 
the use of the non-linear growth to constrain the bias between galaxies and dark matter in 
the same way as it is done in standard case, thus inferring the linear growth factor from 
the normalization of the power spectrum in the linear regime. Although the result in the 
right panel in Fig. 2 may seem a coincidence at first sight, Eq. (3.1) says that the nontrivial 
correction from DGP gravity in square brackets is maximum when 6 = (which gives the 
renormalization of Newton's constant). As 6 increases the correction disappears (since DGP 
becomes Einstein at high-densities), so most of the difference between the two evolutions 
happens in the linear regime, which is encoded in the linear growth factor. 



D. Late-time ISW Effect 



We now consider the late-time integrated Sachs- Wolfe (ISW) effect on the cosmic mi- 
crowave background (CMB) for perturbations with scale r ^ H^^. For this purpose, we 
need to identify the gravitational potentials for linear overdensities as perturbations around 
a homogeneous cosmological background with the line element 

ds^ = [1 + 2$(t, A)] dr^ - d^{r) [1 + 2*(t, A)] [dA^ + X'^dn] . (3.8) 

Here $(t. A) and "^{t, A) are the relevant gravitational potentials and A is a comoving radial 
coordinate. In effect we want to determine $ and ^ given n and a. Unlike the case of 
Einstein's gravity, $ ^ — ^, due to additional contribution of the first term in Eq. (2.4). 
One may perform a coordinate transformation to determine that relationship. We find that, 
assigning r = a(T)A, and 
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background] ) 



(3.9) 
(3.10) 



keeping only the important terms when rH <C 1. But since we are concerned with hnear 
density perturbations, we find from Eqs. (2.23) and (2.24) that the quantity of interest for 
the ISW effect is the time derivative of 



where V is the gradient in comoving spatial coordiantes. This result is identical to the 
four-dimensional Einstein result, the contributions from the brane effects exactly canceUing. 
This result is not entirely surprising since the effect of the brane is the introduction of a 
new gravitational scalar that couples to the trace of the energy-momentum tensor. However, 
the ISW effect has to do with the evolution of photons through a gravitational field (in the 
weak field limit), and photons will not couple to the gravitational scalar (its trace vanishes). 
Thus, the late-time ISW effect for DGP gravity will be identical to that of a dark energy 
cosmology that mimics the DGP cosmic expansion history, Eq. (2.8), at least at scales small 
compared to the horizon. Our approximation does not allow us to address the ISW effect 
at the largest scales (relevant for the CMB at low multipoles), but it is applicable to the 
cross-correlation of the CMB with galaxy surveys. At larger scales, one expects to encounter 
difficulties associated with leakage of gravity off the brane (for order-unity redshifts) and 
other bulk effects [30,29] that we were successfully able to ignore at subhorizon scales. 

Discussion of photon geodesies naturally leads one to ask how lensing may be altered 
due to DGP contributions. For weak lensing by large-scale structure, one is in the weak 
field limit and therefore Eq. (3.11) applies; that is, the weak lensing pattern is identical to 
that for Einstein gravity, apart from the difference in expansion histories and change in the 
force-law. In other words, reconstruction of the dark matter distribution in DGP from weak 
lensing only requires changing the growth rate and the geometrical distances.^ 



^In is a nontrivial result that light deflection by a compact spherical source is identical to that in 
four-dimensional Einstein gravity (even with potentials Eqs. (2.16)-(2.19) substantially differing 
from those of Einstein gravity) through the nonlinear transition between the Einstein phase and 
the weak-brane phase. As such, there remains the possibility that for aspherical lenses that this 
surprising null result does not persist through that transition and that DGP may manifest itself 
through some anomalous lensing feature. 



V^(*-*) 




M2 



(3.11) 
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E. Beyond Isolated Spherical Perturbations 



Since we have derived the growth of spherical isolated perturbations, it is fair to ask how 
well do we expect our results to hold in the realistic case of a superposition of perturbations 
of arbitrary shape. In the linear regime, one expects to recover the same result as here, since 
the linearized equations obey the superposition principle and one may construct arbitrary 
perturbations from a linear superposition of isolated spherical perturbations. In the linear 
regime, DGP gravity reduces to a Brans-Dicke theory with a slowly time-dependent Newton's 
constant, Eq. (3.4). The Newtonian potential is then just a solution to Poisson's equation 
for a given matter distribution source. It would be interesting, however, to corroborate this 
prescription with a fully consistent linear solution for an arbitrary perturbation spectrum 
along the lines of those presented in [30], restricted to scales smaller than the Hubble radius. 

For the non-linear growth, the situation is more complicated since the relation between 
the gravitational potential and the density fluctuation 6 is non-linear and the principle of 
superposition no longer holds. Here, however, one should keep in mind that these non- 
linearities only develop late in the evolution after the universe starts to accelerate, thus the 
corrections to superposition have a small time to act. It is however difficult to say something 
more quantitative at this point. 

IV. OBSERVATIONAL CONSEQUENCES 

What are the imphcations of these results for testing DGP gravity using large-scale struc- 
ture? A clear signature of DGP gravity is the suppressed (compared to ACDM) growth of 
perturbations in the linear regime due to the different expansion history and the addition 
of a repulsive contribution to the force law. However, in order to predict the present nor- 
malization of the power spectrum at large scales, we need to know the normalization of the 
power spectrum at early times from the CMB. A fit of the to pre-WMAP CMB data was 
performed in Ref. [21] using the angular diameter distance for DGP gravity, finding a best 
fit (flat) model with ~ 0.3, with a very similar CMB power spectrum to the standard 
cosmological constant model (with ~ 0.3 and Q° = 0.7) and other parameters kept flxed 
at the same value. Here we use this fact, plus the normalization obtained from the best- 
flt cosmological constant power-law model from WMAP [33] which has basically the same 
(relevant for large-scale structure) parameters as in [21], except for the normalization of the 
primordial fluctuations which has increased compared to pre-WMAP data (see e.g. Fig. 11 
in [34]). The normalization for the cosmological constant scale-invariant model corresponds 
to present rms fluctuations in spheres of 8 Mpc/h, os = 0.9 ± 0.1 (see Table 2 in [33]). 
We assume a flat model, since it was shown in [21] to be consistent as well with the DGP 
angular diameter distance. We ignore the fact that the ISW effect at low multipoles (where 
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FIG. 3. The linear power spectrum normalization, as, for DGP gravity as a function of 
The vertical lines denote the best fit value and 68% confidence level error bars from fitting to 
type-IA supernovae data from [21], = O.lSloQg. The other lines correspond to as as a function 
of obtained by evolving the primordial spectrum as determined by WMAP by the DGP growth 
factor. See text for details. 

the CMB power spectrum has large error bars) can be different in DGP gravity, this has a 
small effect on the overall normalization of the primordial fluctuations that are determined 
by the overall power spectrum. 

Figure 3 shows the present value of erg as a function of for DGP gravity, where we 
assume that the best-flt normalization of the primordial fluctuations stays constant as we 
change Q^, and recompute the transfer function and growth factor as we move away from 
fl^ = 0.3. Since most of the contribution to erg comes from scales r < 100/i/Mpc, we can 
calculate the transfer function using Einstein gravity, since these modes entered the Hubble 
radius at redshifts high enough that they evolve in the standard fashion. The value of erg at 
fl^ = 0.3 is then given by 0.9 times the ratio of the DGP to ACDM growth factors shown in 
the bottom panel of Fig. 1. The error bars in erg reflect the uncertainty in the normalization 
of primordial fluctuations, and we keep them a constant fraction as we vary away from 
0.3. We see in Fig. 3 that for the lower values of fl'^ preferred by fltting the acceleration 
of the universe, the additional suppression of growth plus the change in the shape of the 
density power spectrum drive o"g to a rather small value. This could in part be ameliorated 
by increasing the Hubble constant, but not to the extent needed to keep ag at reasonable 
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values. The vertical lines show the best-fit and la error bars from fitting DGP gravity to 
the supernova data from [1] in [21]. This shows that fitting the acceleration of the universe 
requires approximately ag < 0.7 to la and ag < 0.8 to 2a. 

In order to compare this prediction of ag to observations one must be careful since most 
determinations of as have built in the assumption of Einstein gravity or ACDM models. We 
use galaxy clustering, which in view of the results in Sect. Ill C for higher-order moments, 
should provide a test of galaxy biasing independent of gravity being DGP or Einstein. Recent 
determinations of as from galaxy clustering in the SDSS survey [35] give a^ — 0.89 ± 0.02 
for L* galaxies at an effective redshift of the survey Zg — 0.1. We can convert this value 
to as for dark matter at 2; = as follows. We evolve to 2; = using a conservative growth 
factor, that of DGP for = 0.2. In order to convert from L* galaxies to dark matter, 
we use the results of the bispectrum analysis of the 2dF survey [36] where b — 1.04 ±0.11 
for luminosity L ~ 1.9L*. We then scale to L* galaxies using the empirical relative bias 
relation obtained in [37] that b/b* — 0.85 -|- 0.15(L/L*), which is in very good agreement 
with SDSS (see Fig. 30 in [35]). This implies as = 1.00 ±0.11. Even if we allow for another 
10% systematic uncertainty in this procedure, the preferred value of in DGP gravity that 
fits the supernovae data is about 2a away from that required by the growth of structure at 

An independent way of testing DGP gravity with large-scale structure is to constrain 
the growth of velocity fluctuations through /. This affects the redshift distortions of the 
power spectrum and can be extracted from measurements, though at present the errors 
are somewhat large (see e.g. [35]) for an accurate test, but this should improve soon. The 
interesting feature of this test is that it is independent of the normalization of the primordial 
fluctuations, unlike the as normalization discussed above; thus it will be important to check 
that the same conclusions follow in this case. 

V. CONCLUDING REMARKS 

In this paper we identified how one may test the modifications of the gravitational force 
law expected in Dvali-Gabadadze-Porrati (DGP) gravity at scales of cosmological interest. 
While cosmology is altered when the Hubble scale becomes comparable to today's Hubble 
scale, Hq^, the force law is correspondingly altered at much shorter scales and affects, for 
example, the growth of density perturbations at redshifts of order unity. 

Although the results obtained in this paper are qualitatively comparable to those found 
in [15], they differ in key ways since DGP gravity does not obey Birkhoff's law. The form of 
the gravitational force law between localized mass sources is sensitive to the background cos- 
mological expansion. So, while the results show deviations from Newtonian gravity of order 
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unity at distance scales greater than the scale r* ~ {vg/H^y/^, where Vg is the Schwarzschild 
radius of the mass source, the quantitative details differ. In particular, the suppression of 
the growth of structure is a factor about two larger than in Birkoff-law models with the same 
expansion history. Moreover, DGP gravity deviates significantly from four-dimensional Ein- 
stein gravity through the emergence of an ultra-light graviscalar mode. Because such a mode 
does not couple to photons, the new effects do not manifest themselves in the late-time in- 
tegrated Sachs-Wolfe (ISW) effect, at least at the subhorizon scales we consider. While the 
gravitational potentials are indeed altered by order-unity factors late in the cosmic expansion 
history, they precisely cancel, so that the late-time ISW effect in DGP gravity is identical to 
that for a dark energy theory that mimics the DGP expansion history. The same situation 
applies to weak gravitational lensing. 

We have done a first assessment of the observational viability of DGP gravity to simul- 
taneously explain the acceleration of the universe and the growth of structure. In order to 
improve the comparison against observations a number of issues remain unsolved. First, one 
would like to check that the linear growth factor for sub-horizon scales derived under the 
spherical approximation holds for more general perturbations, as expected by the superpo- 
sition principle in the linear regime. A more non-trivial check would be to generalize this to 
the non-linear case, or at least second-order in perturbation theory, to check the deviations 
from superposition assumed here in the nonlinear case. To do a full comparison of the CMB 
power spectrum against data it remains to solve for the ISW effect at scales comparable 
to the horizon. Our treatment found no difference from general relativity (except from the 
change in the expansion history), but this is only valid at subhorizon scales. 

Nevertheless, the main problem for DGP gravity to simultaneously explain cosmic accel- 
eration and the growth of structure is easy to understand: the expansion history is already 
significantly different from a cosmological constant, corresponding to an effective equation 
of state with Wes = — (1 + flrn)~^- This larger value of w suppresses the growth somewhat 
due to earlier epoch of the onset of acceleration. In addition, the new repulsive contribution 
to the force law suppreses the growth even more, driving erg to a rather low value, in contrast 
with observations. If as error bars shrink the supernovae results continue to be consistent 
with Wes = — 1, this will drive the DGP fit to a yet lower value of and thus a smaller 
value of CTg. For these reasons we expect the tension between explaining acceleration and the 
growth of structure to be robust to a more complete treatment of the comparison of DGP 
gravity against observations. 



17 



ACKNOWLEDGMENTS 



The authors thank C. Deffayet, G. Dvah, G. Gabadadze, A. Gruzinov, M. Takada and 
M. Zaldarriaga for helpful communications and insights. A. L. and G. S. wish to thank 
the CERN Theory Division and the Center for Cosmology and Particle Physics for their 
hospitality. This work is sponsored by DOE Grant DEFG0295ER40898, the CWRU Office 
of the Provost, NASA grant NAG5-12100, NSF grant PHY-0101738, and CERN. G. S. 
thanks Maplesoft for the use of Maple V. 



18 



REFERENCES 



[1] S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 517, 
565 (1999). 

[2] A. G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116, 1009 (1998). 
[3] C. Deffayet, Phys. Lett. B 502, 199 (2001). 

[4] C. Deffayet, G. R. Dvali and G. Gabadadze, Phys. Rev. D 65, 044023 (2002). 

[5] T. Damour, 1. 1. Kogan and A. Papazoglou, Phys. Rev. D 66, 104025 (2002). 

[6] K. Freese and M. Lewis, Phys. Lett. B 540, 1 (2002). 

[7] K. Freese, Nucl. Phys. Proc. SuppL 124, 50 (2003). 

[8] G. Dvah and M. S. Turner, arXiv:astro-ph/0301510. 

[9] Capozziello S., Carloni S., Troisi A., arXiv:astro-ph/0303041 
[10] S. M. Carroll, V. Duvvuri, M. Trodden and M. S. Turner, arXiv:astro-ph/0306438. 
[11] S. Capozziello, V. F. Cardone, S. Carloni, A. Troisi., arXiv:astro-ph/0307018 
[12] S. Nojiri and S. D. Odintsov, Phys. Rev. D 68, 123512 (2003). 
[13] J. Khoury and A. Weltman, arXiv:astro-ph/0309300. 
[14] J. Khoury and A. Weltman, arXiv:astro-ph/0309411. 
[15] A. Lue, R. Scoccimarro and G. Starkman, arXiv:astro-ph/0307034. 
[16] T. Multamaki, E. Gaztanaga and M. Manera, Mon. Not. Roy. Astron. Soc. 344, 761 
(2003). 

[17] G. Dvah, G. Gabadadze and M. Porrati, Phys. Lett. B 485, 208 (2000). 
[18] G. R. Dvah, G. Gabadadze, M. Kolanovic and F. Nitti, Phys. Rev. D 64, 084004 (2001). 
[19] G. R. Dvah, G. Gabadadze, M. Kolanovic and F. Nitti, Phys. Rev. D 65, 024031 (2002). 
[20] P. P. Avehno and C. J. A. Martins, Astrophys. J. 565, 661 (2002); C. Deffayet, 

G. R. Dvah and G. Gabadadze, arXiv:astro-ph/0106449. 
[21] C. Deffayet, S. J. Landau, J. Raux, M. Zaldarriaga and P. Astier, Phys. Rev. D 66, 

024019 (2002). 

[22] J. S. Alcaniz, Phys. Rev. D 65, 123514 (2002); D. Jain, A. Dev and J. S. Alcaniz, Phys. 
Rev. D 66, 083511 (2002); J. S. Alcaniz, D. Jain and A. Dev, Phys. Rev. D 66, 067301 
(2002). 

[23] C. Deffayet, G. R. Dvah, G. Gabadadze and A. 1. Vainshtein, Phys. Rev. D 65, 044026 
(2002). 

[24] A. Lue, Phys. Rev. D 66, 043509 (2002). 

[25] A. Gruzinov, arXiv:astro-ph/0112246. 

[26] M. Porrati, Phys. Lett. B 534, 209 (2002). 

[27] A. Lue and G. Starkman, Phys. Rev. D 67, 064002 (2003). 

[28] G. Dvah, A. Gruzinov and M. Zaldarriaga, Phys. Rev. D 68, 024012 (2003). 

[29] A. Lue, Phys. Rev. D 67, 064004 (2003). 



19 



[30] C. Deffayet, Phys. Rev. D 66, 103504 (2002). 

[31] M. A. Luty, M. Porrati and R. Rattazzi, JHEP 0309, 029 (2003). 

[32] F. Bernardeau, Ap. J. 392, 1 (1992) 

[33] D. N. Spergel et al., Ap. J. S. 148, 175 (2003) 

[34] G. Hinshaw et al, Ap. J. S. 148, 135 (2003) 

[35] M. Tegmark, et al, arXiv:astro-ph/0310725 

[36] L. Verde, et al, MNRAS 335, 432 (2002) 

[37] P. Norberg, MNRAS 328, 64 (2001) 



20 



